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Abstract 

f(R) gravity is one of the simplest generalizations of general relativity, which may explain the 
accelerated cosmic expansion without introducing a cosmological constant. Transformed into the 
Einstein frame, a new scalar degree of freedom appears and it couples with matter fields. In order 
for f(R) theories to pass the local tests of general relativity, it has been known that the chameleon 
mechanism with a so-called thin-shell solution must operate. If the thin-shell constraint is applied 
to a cosmological situation, it has been claimed that the equation-of-state parameter of dark energy 
w must be extremely close to — 1. We argue this is due to the incorrect use of the Poisson equation 
which is valid only in the static case. By solving the correct Klein-Gordon equation perturbatively, 
we show that a thin-shell solution exists even if w deviates appreciably from —1. 
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I. INTRODUCTION 


The cosmic acceleration was discovered from observations of type la supernovae l|, 2], 
The mechanism which causes the acceleration is still not clear and many approaches have 
been tried. The simplest one is the A cold dark matter (ACDM) model, which is based on 
Einstein’s general relativity (GR) and consists of CDM and a cosmological constant A. In 
this model, however, there is a problem in that the parameter must be fine-tuned in order 
to explain the observed energy budget of the Universe. 

Another approach is to modify the theory of gravity from GR, and many models have 
been proposed (for a review, see Ref. 3 ,. One of these is a class of f W theoties M, 
which is the easiest generalization of GR. These models can be recast into the form of GR 
plus a scalar field by a conformal transformation 8 ]. In this Einstein frame, the scalar 


couples with a matter field and matter experiences the fift 
must be small in order to pass the local tests of gravity 


'orce. This new type of force 


9412], To avoid this difficulty, 


it is necessary that the so-called chameleon mechanism operates and a thin-shcll scalar 
configuration exists [ijj]. If an object has a thin shell, the mass of the scalar field gets large 
inside the object and the fifth force is suppressed. As an alternative model for dark energy, 
it is quite important to get a constraint on the effective equation-of-state parameter w from 
the viewpoint of distinguishing models. If the thin-shell constraint is naively applied to 


a cosmological 


situation, it has been claimed that w must be extremely close to — 1 


14] 


(see also Ref. Usj]). However, this constraint is physically unacceptable since a cosmological 
background quantity like w should not be constrained only by local information. In this 


14] 


paper, we perform a more precise analysis and show that the previous constraint of Ref. 
does not apply. 

This paper is organized as follows. In Sec. El we review the basics of f(R) theories. 
The thin-shcll solution, which explains how GR is restored in f(R) gravity in a local scale, 
is also described for a static case. In Sec. m we analyze f(R) cosmology and consider 
the constraint on w. Here we obtain the scalar field configuration with a thin shell in the 
situation where the Universe is dominated by dark energy with w 7 ^ —1. 

Throughout this study, we use natural units with c = h = 1. The reduced Planck mass 
is written as Mpi = ( 87 rG ) -1//2 = re -1 . The sign convention is as follows. 


The metric has signature (—, +, +, +). 
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• The Riemann tensor is defined as R x flU(T = Y x ^ v - T x IXVja + T x au T a firr - Y x aa T a ^ v 

• The Ricci tensor is defined as R /w = R 01 ^. 


II. f(R) GRAVITY AND THE CHAMELEON MECHANISM 

A. /(/?,) theories 

We analyze metric f(R) gravity, whose action is given by replacing the Ricci scalar R in 
the Einstein-Hilbert action by a general f un ction of R: 


S f(R) ~ 


M 


pi 


d A xy/=gf(R) + S m [g T], 


( 1 ) 


where S™ is the action for a matter field T. Note that there exist other formulations of f(R) 


gravity |6|, such as the Palatini formalism 
assumed to be independent variables. 

Variation with respect to g M1/ gives the equation of motion: 


15| . in which the metric and the connection are 


R^uF(R) - - f(R)g = 8vr GTf,,, + V At V„F(i?) - g^DF(R), 


( 2 ) 


where F(R) = df(R)/dR and T jJV is the stress-energy tensor for the matter field. From this 
equation one can see that G e g = G/F{R) acts as the gravitational constant in GR, which 
means Ges depends on both position and time. 

By taking the covariant derivative of Eq. (121) , one can show that 


TV = 0. 


( 3 ) 


So the equation of motion (J2j) is consistent with the conservation of energy-momentum of 
the matter field. 

Now we introduce a scalar field </> by 


F(R) = e 


— p — 2/3</>/Mpi 


( 4 ) 


where /? = l/\/6, and the Einstein-frame metric g )iu is 


g,u = e-W Mp1 ^ 


( 5 ) 
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Then the gravitational part of Eq. (JT]) can be rewritten as 

^erravitv Id X\J (j 


' gravity 

Here V (0) is the potential of 0: 


Mha - ivVv^ - v{4>) 


,,,,, _ RF(R) — f(R) 

VW = ~2 - JW 2 -' 

Note that the right-hand side can be thought of as a function of 0 through Eq. 
In the Einstein frame, the equations of motion are as follows: 


R 




-Rf],w = K V„0V v 0 - K 2 g tw -V a 0 V a 0 + E(0) + nT lwl 


□0 = V(0) - 


M, 


T. 


pi 


( 6 ) 


( 7 ) 


( 8 ) 

( 9 ) 


Here a prime denotes a derivative with respect to 0, T^ v is the stress-energy tensor and 
T = g^'T^ is its trace. T IIV is related to the Jordan-frame quantity by 

fH v = e W/« P1 p^ ( 10 ) 


Since 0 is coupled to the matter held through g jW in the contribution of matter appears 
in Eq. ((2]). Due to this term, T0, is not conserved in the Einstein frame. It is convenient to 
define a conserved quantity, 


_ g -/30/M P1 J'V’v — g3^</>/Mpi 


( 11 ) 


The corresponding densities in the two frames are related by p — e 3 ^^/ Mpi p. For 0 -C Mpi, 
p and p almost coincide with each other. 

Now we have the problem of choosing a physical frame, he., we have to match the Jordan- 
or Einstein-frame quantity with the observed one. In this paper we regard the Jordan frame 
as physical, where the masses of particles do not depend on their position or time. On 
the other hand, it is difficult to solve Eq. (J2]) directly since it is a fourth-order differential 
equation. Thus, we perform calculations in the Einstein frame where the equations of motion 
(jHJ) and (J9J) are both secondorder. 


B. Chameleon mechanism 

In the Newtonian limit, a test mass M experiences a force F given by 0 Q 

F (3 


M 


Mi 


-V0. 


pi 


( 12 ) 
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This fifth force violates the weak equivalence principle and may result in disagreement with 
the local tests of GR. 

In order for f(R) theories to pass the local tests of gravity, the chameleon mechanism 


must work to suppress the fifth force compared to the Newtonian force 


13]. Since we are 


□0 = V'(4>) + Yf = VM), 

M pi 


interested in the present era of accelerated expansion, we neglect radiation and assume the 
Universe is dominated by pressureless matter. The Klein-Gordon (KG) equation (J9]) becomes 

(13) 

where 14 ^( 0 ) is the effective potential, which includes the contribution of the matter field: 

U eff (0) = U(0)+pe^ Mp h (14) 

If the functional form of f(R) is designed to satisfy U'(0) < 0 , U"(0) > 0 and U"'(0) < 0 , 
one can show the following. 

• V e s(4>) has a minimum. 

• The minimum 0 min is a decreasing function of p. 

• The scalar mass at the minimum is an increasing function of p. 

So in a dense region, takes a large value and the range of the fifth force becomes short. 
This is called the chameleon mechanism, and the corresponding scalar is called a chameleon 
field. By virtue of this mechanism, GR is restored on a local scale. 

To see how the fifth force gets small concretely, let us consider a scalar field around 
a uniform spherical object with radius R c , density p c , and mass M c = 47rR]]p c /3. The 
background spacetime has density pb and is assumed to be static. If the object is large 
enough, the value of 0 inside the sphere is given by the minimum of !4ff(0) whose shape_is 
determined by p c . This type of scalar configuration is known as the thin-shell solution 13]. 

The chameleon configuration is given by solving Eq. (fT31) neglecting the time derivatives, 
i.e. the Poisson equation V 2 0 = 15^(0). 0 is assumed to sit at the potential minimum both 
well inside (r < R s < R c ) and far from (r —> oo) the object (0 C and 0&, respectively), and 
the corresponding masses are denoted by m c and mb 21|. The solution is obtained as follows: 

r < R s , 

3 Rl l + 50c, R s <r < R c , ( 15 ) 


50 = < 


50c, 

(3pc 


6M P1 

Ppc 


h + 2 R 

r 


r? 3 

1L C -m b (r-Rc) 

M c th c i 

PI T 


r > R r . 
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Here 5(j) c = (f) c — (j)j, and e t h is a constant which parametrizes the thinness of the shell-like 
region R s < r < R c : 

lA^A I 7? — 7? 

(16) 


_ M P1 |(50 c 

Gh = 


/3 R 2 p c 


/h. 


This means that, for an object to have a thin-shell solution, e t h < 1 is needed. We call this 
inequality the thin-shell constraint. 

Note that e t h can be rewritten as 

/Wcl /M pi 


Hh — 


(17) 


GM C /.R C ’ 

which is considered as the ratio between the fifth force and the Newtonian potential. There¬ 
fore the thin-shell constraint implies that the fifth force is smaller than the Newtonian force. 
From the Solar System tests of the weak equivalence principle using the free-fall acceleration 
of the Earth and the Moon toward the Sun, it is known that eth for the Earth (e t h,®) must 
be smaller than 2.2 x 10 -6 [jj. 


III. EQUATION OF STATE OF DARK ENERGY 


A. f(R) cosmology 


From now on we impose the flat Friedmann-Lemaitre-Robertson-Walker metric in the 
Jordan frame: 

ds 2 = — dt 2 + a(t ) 2 [dr 2 + r 2 (dd 2 + sin 2 6 dip 2 )] . (18) 

Transformed into the Einstein frame, Eq. (jHj) becomes 


H 2 

2d 

a 


1 


Mi 


pi 


H 2 - 


3Mp, L 


W + j], 


W_ 

M 


pi 


+ 2 H6) + 


3Mpj 


A/pi 


V{<f>)F. 


(19) 

( 20 ) 


To interpret these equations in terms of GR, we reassemble Eqs. (HUD an d (EUP as follows: 

H 2 = = HH f! m + n DE ), (21) 

( 22 ) 


-h H 2 = 

a 


3M2jF 0 

Pde 


MpiFo 


where Q m , Ode are the effective density parameters of matter and dark energy, respectively, 

Pm o _ d DE 


Urn - 


3M 2 M| 1 F 0 ’ 


f^DE = 


3// 2 M 2 jF 0 ' 


(23) 
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Equations fl2T]) and fl22|) mean that we regard the present value of the effective gravitational 
constant G e ff as the gravitational constant of GR, and any deviation from GR is considered 
as dark energy: 



(24) 


(25) 


Together with the effective equation of state of dark energy 



(26) 


we get 



(f - i) ■ (27) 


Since time derivatives are of order H , this expression can be estimated as 


|(1 + w)Hde| ~ O 



(28) 


where A0 is the variation of 0 in the last Hubble time. 

Now let us apply the thin-shcll constraint to a cosmological situation, following the steps 
of Ref. jljj]. Here the following assumptions are made. 

(i) The Universe is approximately homogeneous when coarse grained over scales larger 
than some Thom II 1 . 

(ii) Gravity is weak, i.e. <h N <C 1 and v l v l <C 1. 

(iii) The chameleon mechanism works similarly as in the static case. 

Since p c > p b (t ) > pb(to) for a past time t at which z > 1, the following relation holds by 
virtue of the chameleon mechanism: 


0 C < <f) b (t) < 0 b (t O )- 


(29) 


If one identifies 0 &(£q) — 4>bif) as A0 in Eq. ([25]) . it follows that 



( 30 ) 
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On the other hand, the thin-shell constraint says that 


Mi 


o) — (p c ) < 4?n 


(31) 


pi 


where to is the present time. Combining Eqs. (l30l) and (l3Tjh one gets 


1(1 + M7)Ode| < ^N, 


(32) 


where is the Newtonian potential of the celestial object under consideration. Since 
is 10 -6 — 10~ 5 for large clusters and superclusters [lsj], the following constraint is obtained: 


|(1 + w)Ode| < 10 


-4 


(33) 


This can be thought of as a constraint on w itself, since we already know Ode ~ 0.7 19]. 
The constraint fl33l) would imply that f(R ) gravity is indistinguishable from the cosmological 
constant model, as far as the evolution of the homogeneous background is concerned. 

However, this constraint is physically unacceptable because a background cosmological 
quantity like w should not be determined only by local information. This unphysical result 
was derived due to the assumption (iii). Actually, for models which predict m$ 0{H), 
(p does not sit at the minimum of V e g (0) and Eq. (l29]i no longer holds. Furthermore, the 
original exterior solution, which was derived neglecting the time-derivative terms, does not 
satisfy the KG equation f[T3|) . Therefore one must solve the correct held equation for such 
models, which will remedy the constraint of Ref. [bfj]. 


As a specific example, let us adopt Starobinsky’s model 

f(R) = R + \Rt 


20 !]: 



(34) 


with n, A > 0 and R s is of the order of the observed cosmological constant. The parameter 
space can be constrained by both background-level and perturbative-level arguments 20]. 
At the background level, the asymptotic de Sitter solution should be stable. The smallest 
values of A which satisfy the stability condition are 0.95, 0.73, and 0.61 for n — 2, 3, and 4, 
respectively. At the perturbative level, assuming that the ratio of the linear density pertur¬ 
bation 5 /(k)/<5acdm is small, the constraint resembles Fig. 5 of Ref. |20]. In order to keep the 
ratio smaller than 10% at k — 0.174 h Mpc -1 (which is the wave number corresponding to 
(7 g normalization), A should be larger than 8.2, 3.0, and 1.9 for n = 2, 3, and 4 respectively. 
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One can calculate e t h,® in this model and it turns out to be small enough to satisfy the 
Solar System tests of the weak equivalence principle. For example, for n = 2 and A = 1, 

e th,® ~ 10 17 . 

It is known that these models predict |l + iu| ~ 0(0.1) and ~ O(H). So if the solution 
for the chameleon held with a thin shell is obtained in this framework, it can be thought of 
as a counterexample of the previous work {l^ . 

B. Chameleon configuration in the accelerating Universe 

From now on we calculate a configuration for the chameleon field around a spherical 
object in a case where the deviation of w from —1 is fairly large. The following assumptions 
are made. 

• The background spacetime evolves as in the wCDM model (CDM+dark energy with 
constant w). 

• The object is decoupled from cosmic expansion and its radius R c is constant in physical 
coordinates, i.e., its surface is at ar = R c . 

• The interior (ar < R c ) solution 8(f)- m is given by Eq. (fT5jh in which r is replaced by ar. 

• The exterior (ar > R c ) solution 8(p out is smoothly connected to h0 in . 

• 5(j) approaches zero at infinity. 

Here 8(p denotes the deviation of </> from the cosmological background value 4>b{t), which 
does not necessarily correspond to the minimum of V e s((p). 

First we consider the de Sitter (w = —1) case, for which we can obtain the solution 
analytically, followed by the w ^ — 1 case. In the following we use conformal time rj as the 
time variable. 

1. w = — 1 case 

The KG equation outside of the object becomes 

— 84>" — 2 R8(f)' + V 2 8(j) — m 2 a 2 8(j) = 0, (35) 
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where % = a'/a, a prime denotes a derivative with respect to 77, and mb is the background 
mass of ep. Since ep = —(Mpi/ 2 / 3 ) In f'(R) = const in de Sitter spacetime, we neglect the 77 
dependence of the background ep and mb- Using a = —( Hr;) _1 for w = —1, Eq. (1H51) can be 
written as follows: 

9 jy) 2 

- Set)" + -Set)' + V 2 Sej) - -j^Set) = 0. (36) 

77 H z r] z 

The boundary conditions are 

• <50out Sep in at (XT Rci 

• d r Sep out =d r Sepi n at ar = R CJ 

• Sep out —» 0 as ar —> 00 . 

Let us describe the solution Sep in terms of conformal time and the physical distance 
instead of the comoving distance, be., in the form 


S<p(rj,r) = 


(37) 


where u denotes the physical distance normalized by H h 

r 


u = Har = T~Lr = — - 


77 


(38) 


Then Eq. 


becomes 


02 2 o 1 

C'nA C't/C'uA T (Xrj^P 9 


(u 2 — 1 )<9„a + — — 'j d u (fi + ^ 


f m b y 


A 


= 0. (39) 


The expression inside the square brackets of Eq. (1391) is written only with w. Also note that 
the boundary conditions for ip are 

® Aout Ain ^ dd R c: 

* ^wAout S) u (pin &t U HR Cl 

• Aout —* 0 as u —> 00 , 

and are expressed only with u. So we investigate a solution that is independent of 77 , be., 
tp = ip(u ). In such a case, Eq. (1391) can be rewritten as the following ordinary differential 
equation: 

+ 4u i - 2 , a + = 0 . (40) 


du 2 u(u 2 — 1 ) du 


u 2 — 1 
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The independent solutions are obtained in terms of the hypergeometric function 2 T\ as 
follows: 




3 + 2ia 3 — 2 ia 3 2 ' 

4 ’ 4 '2' U 

f 1 + 2 ia 1 — 2 ia 1 
\ 4 ’ 4 ’ 2 ’ U 


(41) 

(42) 


where 



(43) 


Both (pa\u) and <pa\u) diverge at u — 1. Fortunately, we can construct a solution g a (u) 
which is finite at u = 1 by taking their linear combination: 


9a(u) = ip£\u)-2 


p ^ 3+2 ia ^ p ^ 3—2 ia ^ 
p ^ 1+2 ia ^ p ^ 1—2 ia 



(44) 


Note that this function is defined for u > 0 and takes a real value. The general form of 
g a {u) is shown in Fig. |T| 



FIG. 1: The exterior solution g a {u) for a = 1. The dashed line is a plot of 1/u. 


The exterior solution is proportional to g a {T-Lr ), and the factor of proportionality is de¬ 
termined from the boundary conditions at the surface of the object ar = R C1 or u = HR C . 
Using the fact that g a (u ) can be approximated 22] as 1/u for u ~ HR C 1, we get the 
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following configuration for 5(j): 


( 


8(/>c, 

Ppc 


(ar) 2 + 2— — 3 R 2 + S(j) c , R s < ar < R c , 


ar < R s , 


^ - 6 Afpi 


ar 

>2 


(45) 


Mpi 


e th HR c g a (Hr ), ar > R c . 


ar > R c . 


Here e t h is given by Eq. (1XH[) . Thus we managed to find the exterior solution that satisfies 
the boundary conditions. 

2. w A — 1 case 

Now we move on to the case where w is a constant but not —1. We write w = — 1 + e 
and consider up to first-order terms in e. In this case, R is calculated as 



(46) 


Therefore the equation of motion gets slightly modified from Eq. (j36]h and so does the 
solution. We assume the following form for the exterior solution: 


H = - l3pcRc e th HR c [g a (Rr) + eA(r], r)]. 
M pi 


(47) 


This expansion is valid if A < 0(1) at the horizon scale, since g a (Rr) ~ 0(1) there (Fig. [T]). 
The equation for the perturbative part A can be written as 


e[A" + 2 HA! - V 2 A + m 2 b A} = 4 [-(2C 0 - 3 )Urg' a (Ur) - 2C 4> g a {Hr )}, (48) 

rj 


where C^ characterizes the rate of change of (J) b , 


Here C$ ~ 3n for Starobinsky’s model. 

Again we try to find a solution in the form 



(49) 


A(rj,r) = B(u), u ——. 

V 


r 


(50) 


Then Eq. (1481) is rewritten as an ordinary differential equation for B(u ): 



( 51 ) 
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where the source term is given by 


• / S _ ( 2C V - 3)ug' a (u) + 2CVflq(u) 

M 2 — 1 

Since we know the homogeneous solutions for Eq. f!51[) . the inhomogeneous solutions can be 
obtained by the method of variation of parameters. We choose 


Bi(u) = <p£\u), 

(53) 

B 2 (u) = g a (u) 

(54) 


as a basis for the vector space spanned by the homogeneous solutions. Using these, the 
inhomogeneous solution for Eq. d5Th can be written as follows: 


B{u) = CiBiiu) + C 2 B 2 (u) + B s (u), (55) 

B s (u) = —Bi(u) dx^^j{x) + B 2 (u) dx ^^j{x) (56) 

where Ci,C 2 are constants and W[u) is the Wronskian of B\{u), B 2 {u). The coefficient C\ 
is fixed by assuming B(u) does not diverge at u — 1: 

Ci = l dx Wi jix) ■ (57) 

C 2 is determined by requiring B{u) =0 at the surface of the object u = HR C \ 


C 2 = 


Bi{HRJ I 

B 2 (HR c ) Jhx~~W(x) 


dx?Mj(x) - 


-HR C 


dx Bl ^ j(x) 
W(x) J[ } 


(58) 


Combining these, we get the solution for the perturbative part: 


B{u) = —Bi{u) 


dx WWj j( - x)+B2 ^ 


r dx^iux, 

JHR W(X) J(X> 


Bi(HR c ) I 

B 2 (HR c ) J hr ~~W(x) 


dx?Mj(x) 


(59) 

The form of B(u) for various parameters of Starobinsky’s model is shown in Fig. [21 Here we 
set HR C = 1CU 3 , which corresponds to the ratio of the typical scale of a galaxy cluster to 
the present Hubble radius 23]. Even in the n — 2, A = 1 case, where the deviation from GR 
is the largest among the viable models, the magnitude of B(u) is < 1.3. Thus we managed 
to obtain a consistent perturbative solution in the case of w ^ — 1. 
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FIG. 2: The perturbative part B(u) for various parameters of Starobinsky’s model. 


IV. CONCLUSION 


Transformed into the Einstein frame, f(R ) gravity predicts a fifth force due to a scalar 
which is nonminimally coupled to matter. The fifth force is small if an object has a thin- 
shell configuration of the scalar field, and this enables f(R ) theories to pass local tests of 
gravity. If the thin-shell constraint is naively applied to cosmological scales, it leads to an 
extremely small |1 + w\. This result was derived by solving the Poisson equation for the 
scalar held, which is inappropriate because the evolution of the scalar held is determined 
by the Klein-Gordon equation. By solving the correct held equation, we have shown that a 
consistent solution exists even if w deviates appreciably from—1 as long as other viability 
conditions are satished. 
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